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Abstract 

We study the associated production of bottomonia and Higgs bosons at hadron 
colliders within the factorization formalism of nonrelativistic quantum chromody- 
namics providing all contributing partonic cross sections in analytic form. While 
such processes tend to be suppressed in the standard model, they may have interest- 
ing cross sections in its minimal super symmetric extension, especially at large values 
of tan /3, where the bottom Yukawa couplings are enhanced. We present numerical 
results for the processes involving the lighter CP-even h boson and the CP-odd 
A boson appropriate for the Fermilab Tevatron and the CERN LHC. 
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1 Introduction 



Since the discovery of the J/i[) meson in 1974, heavy quarkonium has provided a useful 
laboratory for quantitative tests of quantum chromodynamics (QCD) and, in particular, 
of the interplay of perturbative and nonperturbative phenomena. The factorization for- 
malism of nonrelativistic QCD (NRQCD) [1] provides a rigorous theoretical framework 
for the description of heavy-quarkonium production and decay. This formalism implies 
a separation of short- distance coefficients, which can be calculated perturbatively as ex- 
pansions in the strong-coupling constant a s , from long-distance matrix elements (MEs), 
which must be extracted from experiment. The relative importance of the latter can be 
estimated by means of velocity scaling rules; i.e., the MEs are predicted to scale with a 
definite power of the heavy-quark (Q) velocity v in the limit t) < 1. In this way, the 
theoretical predictions are organized as double expansions in a s and v. A crucial feature 
of this formalism is that it takes into account the complete structure of the QQ Fock 
space, which is spanned by the states n = with definite spin S, orbital angular 

momentum L, total angular momentum J, and color multiplicity a = 1,8. In particular, 
this formalism predicts the existence of color-octet (CO) processes in nature. This means 
that QQ pairs are produced at short distances in CO states and subsequently evolve into 
physical, color-singlet (CS) quarkonia by the nonperturbative emission of soft gluons. In 
the limit v — > 0, the traditional CS model (CSM) [2] is recovered. The greatest triumph 
of this formalism was that it was able to correctly describe [3-5] the cross section of inclu- 
sive charmonium hadropro duct ion measured in pp collisions at the Fermilab Tevatron [6] , 
which had turned out to be more than one order of magnitude in excess of the theoretical 
prediction based on the CSM. 

Apart from this phenomenological drawback, the CSM also suffers from severe concep- 
tual problems indicating that it is incomplete. These include the presence of logarithmic 
infrared divergences in the 0(a s ) corrections to P-wave decays to light hadrons and in 
the relativistic corrections to S'-wave annihilation [7], and the lack of a general argument 
for its validity in higher orders of perturbation theory. While the /^-factorization [8] and 
hard-comover-scattering [9] approaches manage to bring the CSM prediction much closer 
to the Tevatron data, they do not cure the conceptual defects of the CSM. The color evap- 
oration model [10], which is intuitive and useful for qualitative studies, also significantly 
improves the description of the Tevatron data as compared to the CSM [11]. However, 
it does not account for the process-specific weights of the CS and CO contributions, but 
rather assumes a fixed ratio of 1 : 7. In this sense, a coequal alternative to the NRQCD 
factorization formalism is presently not available. 

In order to convincingly establish the phenomenological significance of the CO pro- 
cesses, it is indispensable to identify them in other kinds of high-energy experiments as 
well. Studies of charmonium production in ep photoproduction, ep and uN deep-inelastic 
scattering (DIS), e + e~ annihilation, 77 collisions, and 6-hadron decays may be found in 
the literature; see Refs. [12,13] and references cited therein. Furthermore, the polariza- 
tion of ip' mesons produced directly [14] and of J/ip mesons produced promptly [15,16], 
i.e., either directly or via the feed-down from heavier charmonia, which also provides a 
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sensitive probe of CO processes, was investigated. Until recently, none of these studies 
was able to prove or disprove the NRQCD factorization hypothesis. However, HI data of 
ep -> e J/ip + X in DIS at HERA [17] and DELPHI data of 77 -> J/ip + X at LEP2 [18] 
now provide first independent evidence for it [19,20]. 

Recently, we studied the associated production of charmonia and electroweak bosons 
at present and future colliders, including the Tevatron (Run II), the CERN Large Hadron 
Collider (LHC), the DESY TeV- Energy Superconducting Linear Accelerator (TESLA), 
and the DESY electron-proton supercollider THERA, which uses the lepton beam of 
TESLA and the proton beam of the DESY Hadron-Electron Ring Accelerator (HERA) 
[21]. Our assessment of the observability of the various processes in the various experi- 
ments can be summarized as follows: the processes with a photon will abundantly take 
place in all considered experiments; the processes with a Z or W boson will produce con- 
siderable yields at the hadron colliders, namely, several hundred (ten thousand) events per 
year at the Tevatron (LHC), while they significantly fall short of the one-event-per-year 
mark at TESLA and THERA; the processes with a standard model (SM) Higgs boson 
(H) are predicted to be too rare to be observable in any of the considered experiments. 
As for J/ip + W associated production at the Tevatron, the conclusions that had been 
reached in Ref. [22] concur with ours [21]. The associated production of bottomonia and 
Z or W bosons at the Tevatron and the LHC was investigated in Ref. [23], with the 
conclusion that such events represent a challenge at the Tevatron, while they should cer- 
tainly be observable at the LHC. The CDF Collaboration [24] did not find evidence for 
T(1S) + Z or T(IS') + W ± signals in Run I at the Tevatron, but established upper bounds 
on the respective production cross sections, which are well compatible with the theoret- 
ical predictions of Ref. [23]. However, with the expected Run II increase of integrated 
luminosity and the upgrades of the CDF detector, they expect to achieve a sensitivity 
sufficient to observe these signals [24]. The cross sections of all partonic processes con- 
tributing to the associated production of heavy quarkonia, with 2S+1 Lj = ^o, 3 Si, ^1, 3 Pj 
(J = 0,1,2), and photons, Z bosons, and W bosons in photon-photon, photon-hadron, 
and hadron-hadron collisions may be found in Ref. [21]. 

In this paper, we take the next step and consider the production of bottomonia in 
association with Higgs bosons concentrating on hadron-hadron collisions. Guided by our 
previous findings in connection with charmonium [21], we expect such signals to be below 
the edge of observability in two-photon collisions at TESLA and photon-hadron collisions 
at THERA. Therefore, we do not include the latter production mechanisms in our present 
phenomenological analysis. However, the relevant partonic cross sections may be obtained 
from our analytic results by adjusting overall color factors and coupling constants, as is 
explained in the Appendix. 

In the SM, the total cross section of pp — > J/ip + 7i, with Higgs-boson mass = 
115 GeV, at the LHC was predicted to be 2.0 x 10 -2 fb, which corresponds to 4 signal 
events per year [21]. Taking into account the branching fractions of the J/vp and 7i decays 
to the detected particles and the acceptance cuts required for background suppression, it 
becomes apparent that the production rate would need to be increased by a few orders of 
magnitude before a sound signal could be observed. Since the SM Yukawa couplings are 
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proportional to the quark masses, passing from charmonium to bottomonium yields an 
enhancement by a factor of m 2 /m 2 ~ 10. Moreover, the MEs of bottomonium [13,25] are 
typically one order of magnitude larger than their counterparts for charmonium [15], with 
the exception of (0 T< ™ 5 )[ 1 So 8 ' ) ]) and (C^'^ppj 8 ^), which tend to be comparable to their 
counterparts for ip(nS). On the other hand, the dynamical quark- mass dependence acts 
in the opposite direction. In fact, in the case of the total cross section of pp — > J/ip + H, 
with rrifi = 115 GeV at the LHC, replacing m c by m b while keeping the Yukawa coupling 
and the MEs fixed leads to a decrease by a bit more than one order of magnitude, so that 
the Yukawa coupling enhancement is slightly overtaken. Detailed calculation reveals that 
the total cross section of pp — > T(IS') + 7i at the LHC exceeds the one of pp — > J/ip + H 
by a factor of 3.6, the result being 7.3 x 1CT 2 fb, which corresponds to 15 signal events 
per year. 

The situation should be considerably more promising in the context of supersymmetric 
theories, where the Yukawa couplings depend on additional input parameters and may 
be significantly enhanced relative to the reference values in the SM if nature has chosen 
these parameters accordingly. In the following, we concentrate on the most popular of such 
theories, the minimal supersymmetric extension of the SM (MSSM). The Higgs sector of 
the MSSM consists of a two-Higgs-doublet model of type II and accommodates a quintet 
of physical Higgs bosons: the neutral CP-even h and H bosons, the neutral CP-odd A 
boson, and the charged PT ± -boson pair. At the tree level, the MSSM Higgs sector has 
two free parameters, which are usually taken to be the mass tjia of the A boson and 
the ratio tan/5 = V2/V1 of the vacuum expectation values of the two Higgs doublets. The 
masses rrih, fiiff, and m#± of the h, H, and bosons, respectively, and the mixing angle 
a that rotates the weak CP-even Higgs eigenstates into the mass eigenstates h and H 
are then functions of and tan (3. As is well known, these functions receive significant 
electroweak radiative corrections, whose leading terms are of O (GFmf/m 2 z ), where Gp is 
Fermi's constant, m t is the top-quark mass, and m z is the Z-boson mass, and must not be 
neglected; for a review, see Ref. [26]. For the sake of our exploratory study, it is sufficient 
to make a few simplifying assumptions in the treatment of these corrections. Specifically, 
we neglect subleading effects due to nonzero values of the Higgs-higgsino mass parameter 
fi and the trilinear Higgs-sfermion couplings A t and A b . The leading radiative corrections 
can then be parameterized by the quantity 




(1) 



where ms is the common squark mass, and we have 



2 x 2 1 2 1 
m h = o m A + m z + e 



— y { m \ + m 2 z + e ) 2 — 4m^m| cos 2 (2/3) — 4e (m 2 A sin 2 f3 + m| cos 2 (3^j , (2) 

2 2,2 2 1 fo\ 

m H =m A + m z -m h + e, (3) 



4 



m H ± =m A + m w , (4) 

tan(2a) -(mi-m|)cos(2/3)+e' (5) 

where mw is the H^-boson mass. Solving Eq. (5) for a, one needs to select the branch 
with — 7r/2 < a < 0. When m ft and tan/? are chosen to be the input parameters, then 
Eq. (2) can be solved for tua, which leads to 

2 m^(m| -m 2 h ) + e {m\ - m| cos 2 /?) 
^ m|cos 2 (2/3)-m 2 +esin 2 /3 ' {> 

and a can be obtained from Eq. (5) as before. From Eq. (6), we glean that the upper 
bound on m\ is shifted from its tree-level value, m\ cos 2 (2/3), to 

m 2 h < m 2 z cos 2 (2/3) + e sin 2 (3. (7) 

Combining Eqs. (2) and (3) and noticing that the radicant in Eq. (2) can be written as 
[m\ - m| cos(4/3) + ecos(2/3)f + [m| sin(4/3) - esin(2/3)] 2 , we learn that 

m 2 H > ml + 



mf sin(4/3) -e sin (2/3) . (8) 

Furthermore, inserting Eq. (7) into Eq. (3), we find that 

m 2 H >m 2 A + m| sin 2 (2(3) + e cos 2 (3. (9) 

In other words, H is always heavier than h and A. 

If we assume only the first four quark flavors to be active in the colliding hadrons, 
then the associated production of bottomonium and a charged Higgs boson is suppressed 
by the smallness of the duH~ and scH~ couplings. If we also take bottom to be an active 
flavor, then the final state must contain an additional top quark, W boson, or charged 
Higgs boson, which substantially squeeze the available phase space and lead to distinctive 
signals. For these reasons, we only consider the neutral MSSM Higgs bosons $ = h, H, A 
in the following. The gg<3> couplings emerge by scaling the SM Yukawa couplings with 
factors g®; i.e., 

yf = 2 1 /*G 1 J 2 m q gl (10) 

The values of g® are specified in Table 1. From Table 1, we observe that the values of gf 
rise linearly with tan/3 if tan/3 ^> 1. In order to obtain sufficiently large cross sections, 
we are particularly interested in light Higgs bosons with strong couplings to the b quark. 
We thus focus our attention on the h and A bosons, and on the large-tan (3 regime. 

This paper is organized as follows. In Sec. 2, we present our analytic results and 
explain how to evaluate the cross sections of the associated production of bottomonia 
and Higgs bosons in hadronic collisions. The contributing partonic cross sections are 
collected in the Appendix. In Sec. 3, we present our numerical results. Our conclusions 
are summarized in Sec. 4. 
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Table 1: Values of in Eq. (10) for generic up and down quarks, U and D, respectively. 



1 


9 h a 




3« A 


u 


cos a/ sin f3 


sin a/ sin j3 


cot /3 


D 


— sin a/ cos (5 


cos a/ cos /3 


tan/3 



2 Analytic results 



In this section, we present our analytic results for the cross sections of the reactions 
AB — > CD, where A and 5 are the incoming hadrons, C is a heavy-quarkonium state, 
with 2S+1 Lj = X 3 ^, 1 / 3 !, 3 / 3 / (J = 0,1,2), and D = h, A. The results for D = H, H 
may be obtained from those for D = h by accordingly replacing yq and m^. We also list 
the formulas for the cases where A and/or B are real or quasireal photons, appropriate 
for two-photon and photon-hadron collisions, respectively. 

We work in the fixed-flavor-number scheme; i.e., A and B contain n/ = 4 active quark 
flavors q — u,d,s,c if Q — b and rif — 3 active quark flavors q — u, d, s if Q — c. As 
required by parton-model kinematics, we treat the quark flavors q as massless. The quark 
flavor Q, with mass m,Q, only appears in the final state. The bb Fock states contributing 
at LO in v are n = ^ , X S^ M? M %) if H = r} b (nS); n = 3 S[ 1] , ^ 8) , 3 ^ 8) , 3 Pj 8) if 
H = T(nS); n = iP^, % s) if H = h b (nP); and n = 3 Pf ] , 3 S{ 8) if H = Xbj(nP), where 
J = 0, 1, 2. Their MEs satisfy the multiplicity relations 



Qt{nS) 
QXb.j(nP) 
QXbj{nP) 



3p(8) 
3p(l) 
3^(8) 



(2J+ 1) (C T (" 5 ) 
(2J+l)(c too(riP) 



3p(8) 
-t n 



3 D (1) 



(2J+1)(C» 0( ™ P) fe{ 8) 



(11) 



which follow to LO in v from heavy-quark spin symmetry. The CS MEs of T(nS) and 
Xw(nP) are related to the respective radial wave functions at the origin as 



QXbo(nP) 



3 sf> 



3d(1) 



: ^l^5(0)| 2 , 



;i2) 



respectively. In our numerical analysis, we only include the T(lS') and x &J (lP) mesons. 
For completeness and future use, we also list formulas for all the other bottomonia. The 
assignments for the various charmonia are analogous. 

We now turn to the partonic subprocesses ab — > QQ[n]D. The differential cross 
section of such a process is calculated from the pertaining transition-matrix element T 
as da/dt = |T| 2 /(167rs 2 ), where the average is over the spin and color degrees of freedom 
of a and b and the spin of D is summed over. We apply the covariant-projector method 
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of Ref. [27] to implement the QQ Fock states n according to the NRQCD factorization 
formalism. 

The following partonic subprocesses contribute to LO in a s and v: 

gg^QQi^D, (13) 

gg -> QQ[^ 8) ]D, (14) 

gg^QQ^p, (15) 

ig ^QQ[^)]D, (16) 

QQk (1) ]A (17) 
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where q— 1 Sq, 3 S 1: 1 P 1 , 3 Pj with J = 0, 1, 2. For the reason explained above, g = c must not 
be included if Q — c. The processes — > QQ^WjD and 77 — > <5Q[^®]-D are forbidden 
by color conservation. Furthermore, the processes — > Q<5fc^]-D are prohibited because 
the Q-quark line is connected with the g-quark line by one gluon, which transmits color to 
the QQ pair. Finally, due to charge-conjugation invariance, the processes gg — ► QQ[q^]D, 
ig —> QQ[^]D, and 77 — > QQ[q^]D are forbidden for q = 3 S 1 , 1 P 1 , and the processes 
qq -> QQ[^ 8 )]L> are forbidden for ? = 3 Pj. 

The differential cross sections da/dt of processes (13)— (17) are listed in the Appendix. 



3p(8) 



and 



We combine the results proportional to the CO MEs (^O r ^ nS ^ 
(QXbj(nP) 3gi(8) ^ f Qr g xec [ value of the principal quantum number n, exploiting the mul- 
tiplicity relations of Eq. (11). The results for processes (16) and (17) are obtained from 
those for processes (13) by adjusting overall color factors and coupling constants as ex- 
plained in the Appendix. Similar relations also exist between processes (13) and (14) for 

The evaluation of the hadronic cross sections and their distributions proceeds as ex- 
plained in Ref. [21]. Specifically, we consider the distributions in the transverse momentum 
Pt common to C and D, the rapidities yc and yo of C and D, respectively, and the CD 
invariant mass m co . 



3 Numerical results 

We are now in a position to explore the phenomenological consequences of our calculations. 
We focus our attention on the cases C = Y(lS),XbjO-P)- These bottomonia can be 
efficiently identified experimentally, and their MEs are relatively well constrained [13,25]. 
The predicted cross-section distributions for the T(2S) and T(3S f ) mesons are similar to 
those for the T(1S) meson, but their normalizations are somewhat suppressed due to 
smaller MEs [13,25]. The MEs of the Xb.j(2P) mesons are very similar to those of the 
Xbj{}-P) mesons [13,25], and so are their predicted cross-section distributions. The rjb(lS) 
meson needs experimental confirmation, while no events for r]b(nS) mesons with n > 1 
and hb(nS) mesons have been seen so far [28]. 

We first describe our theoretical input and the kinematic conditions. We use m& = 
4.88 GeV [13,25], m t = 174.3 GeV, m z = 91.1876 GeV, G F = 1.16639 x 1(T 5 GeV~ 2 , 
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and the LO formula for ai nf \fj, r ) [28] with rif — 4 and asymptotic scale parameter 
Aq^ d = 192 MeV [29]. As for the proton PDFs, we use the LO set from the CTEQ 
Collaboration (CTEQ5L) [29]. Unless otherwise specified, we choose the renormalization 



and factorization scales to be /i r = /if = ym^m^, where m% = ym^ +Pt is the trans- 
verse mass of C and similarly for D. We evaluate (o T ^ [ 3 S[ 1] ]) and (p Xb ^ lp "> ^Pq ] \) 

from Eq. (12) with the values of |i?is(0)| 2 and |-R' 1P (0)| 2 determined in Ref. [25] using the 
Q CD- motivated potential by Buchmiiller and Tye [30]. As for the T(IS') and Xbj^-P) CO 
MEs, we adopt as our default set the one determined in Ref. [13] through a fit to recent 
CDF data [31] using the same proton PDFs and CS MEs as we do. Our MSSM benchmark 
scenario is defined by the parameter set tan/3 = 50, m h = 100 GeV or rriA = 100 GeV, 
and ms = 1 TeV, and is well inside the currently allowed region of the MSSM parameter 
space [28]. We then vary tan/3, m^, and itla one at a time in the ranges 2 < tan/3 < 60, 
90 < m h < 128 GeV, and 90 < m A < 500 GeV, respectively. If D — h, then we take 
m h as an input parameter and evaluate uia from Eq. (6). The hadronic center-of-mass 
energy is y/S = 2 TeV in Run II at the Tevatron and = 14 TeV at the LHC. 

Since our study is at an exploratory level, we refrain from presenting a full-fledged 
quantitative estimate of the theoretical uncertainties in our predictions. Experience from 
previous analyses of charmonium production within the NRQCD factorization formalism 
[19,20,32] leads us to expect relative errors of the order of ±50%. This error estimate 
should be on the conservative side, since the convergence property of the NRQCD pertur- 
bative expansion is expected to be considerably improved as one passes from charmonium 
to bottomonium, thanks to the reduction in size of a s and v. To be specific, we assess 
the theoretical uncertainties arising from the lack of knowledge of the precise values of 
the bottomonium MEs and from the freedom in the choice of the renormalization and 
factorization scales. 

We are now in a position to present our numerical results. Figures 1 and 2 are 
devoted to pp — > CD + X at the LHC and to pp — > CD + X in Run II at the Tevatron, 
respectively. In each figure, part (a) gives the pt distributions da/dpr (upper panel) 
and the mcD distributions da/dmco (lower panel), part (b) the yc distributions da/dyc 
(upper panel) and the yo distributions da/dyo (lower panel), and part (c) the total 
cross sections a as functions of tan/3 (upper panel) and mo (lower panel). In each 
part, the left column refers to D = h and the right one to D — A. In each frame, we 
separately consider C = T(1S), x c j(1P), both in NRQCD and the CSM. It is summed 
over C = Xco(l-P), Xci(l-P), Xc2^-P)- In the following, n — 1 is implied whenever the label 
(nL) is omitted. 

As is evident from Sec. 2, T + h and T + A associated production proceeds through 
process (14) with q = ^o, 3 S±, 3 Pj and process (15) with = 3 Si (see solid lines in Figs. 1 
and 2) and is forbidden to LO in the CSM. On the other hand, Xbj + h and Xbj + A 
associated production is already possible in the CSM, through process (13) with q = 3 Pj 
(see dotted lines in Figs. 1 and 2); in NRQCD, also the CO processes (14) and (15) with 
q = 3 Si contribute. 

In the study of the associated production of heavy quarkonia and electroweak gauge 
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bosons [21], dominant contributions at pr 3> rric were found to generally arise from 
so-called fragmentation-prone partonic subprocesses, which contain a gluon with small 
virtuality, q 2 = itiq, that splits into a QQ pair in the Fock state n = Such processes 

are absent here because the Higgs bosons are always radiated from the outgoing Q or Q 
quark lines before the latter form an asymptotic Fock state n. 

We now turn to the predictions for the LHC. From Figs. l(a)-(c), we observe that 
there is hardly any difference between the various distributions for D = h and their 
counterparts for D = A, except for a trivial difference in the mh and tua dependences, 
which is due to the upper bound on m h , at m h 128 GeV. In each frame, the three curves 
exhibit similar shapes (except for the pr distributions), but significantly differ in their 
overall normalizations. Within NRQCD, the results for C = Xbj typically exceed those 
for C = T by a factor of 2. As for C = Xbj, the CSM predictions fall short of the NRQCD 
ones by a factor of 3 to 4. From Fig. 1(a), we read off that the px and uicd distributions 
peak at pr ~ 6 GeV and m C D ~ 116 GeV, respectively. Beyond their peaks, the CSM 
predictions for C = Xbj fall off more rapidly with increasing value of p T than the NRQCD 
ones. The yc and yo distributions in Fig. 1(b) are symmetric about the origin, reflecting 
the symmetry of the experimental set-up, and their maxima are rather broad. Figure 1(c) 
displays the quadratic tan/3 dependence discussed in Sec. 1, which is exact for D = A, 
but only approximate for D = h. For m h <, 120 GeV, the m h dependence is approximately 
exponential, following the heuristic law a(pp — > Ch + X) oc exp[— A(m h /100 GeV)]. On 
the other hand, the dependence is approximately power-like over the range of values 
considered, a heuristic law being a(pp — > CA + X) oc m~ A h . 

We now move on to the predictions for Run II at the Tevatron. Comparing Figs. 2(a)- 
(c) with Figs. l(a)-(c), we observe that all distributions are approximately scaled down by 
a factor of 50 as we pass from the LHC to the Tevatron. Apart from that, their qualitative 
features essentially remain the same. 



Table 2: Total cross sections a in fb of pp — > CD + X at the LHC and of pp — > CD + X 
in Run II at the Tevatron, where C = T(1S), XbjO-P) an d D = h,A, for tan [3 = 50 and 
tjid = 100 GeV. The values enclosed in parentheses refer to the CSM. 



C + D 


LHC Tevatron 


T + h 

XbJ + h 
T + A 
Xb.j + A 


3.3 x 10 2 6.0 

7.4 x 10 2 (1.9 x 10 2 ) 1.3 x 10 1 (3.6) 

3.4 xlO 2 6.1 

7.4 x 10 2 (1.9 x 10 2 ) 1.3 x 10 1 (3.6) 



We now assess the observability of the various processes at the LHC and in Run II 
the Tevatron. To this end, we list their total cross sections for our MSSM benchmark 
scenario in Table 2. These values can be converted into annual yields by recalling from 
Table I of Ref. [21] that a cross section of 0.005 fb (0.25 fb) corresponds to one event 
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per year at the LHC (Tevatron). At hadron colliders, the T(IS') mesons are easily de- 
tected through their decays to e + e~ and fi + fi~ pairs, with a combined branching frac- 
tion of B(T(1S) -> l + l~) = 4.86 ± 0.13 [28]. The XbjO-P) mesons radiatively decay to 
T(IS') mesons, the individual branching fractions being B(xbo0-P) — y ^(1>S , )7) < 6%, 
B( X bi(lP) -> T(1S) 7 ) = (35 ±8)%, and B( Xb2 (lP) -> T(1S) 7 ) = (22 ± 4)% [28]. In 
the high-tan /3 regime, light h and A bosons predominantly decay to bb and t + t~ pairs, 
with branching fractions of about 90% and 8%, respectively [28]. The b hadrons can be 
detected by looking for displaced decay vertices with dedicated vertex detectors, even 
at the LHC. The corresponding numbers of signal events per year estimated assuming 
detection efficiencies of 100% are listed in Table 3. 



Table 3: Annual numbers of events of pp — > CD + X at the LHC and of pp — > CD + X 
in Run II at the Tevatron, where C = T(1S), XbjO-P) an d D = h,A, for tan/? = 50 
and mp = 100 GeV. The final-state particles are assumed to be detected via their decays 
Xbj(lP) -> T(1S)7, T (IS) -> e + e-,n + n~, and h,A^ bb, t + t~ with efficiencies of 100%. 
The values enclosed in parentheses refer to the CSM. 



C + D 


LHC 




Tevatron 


T + h 


3.1 x 10 3 




1.1 


Xbj + h 


1.7 x 10 3 (4.4 x 


10 2 ) 


6.2 x lO" 1 (1.6 x lO" 1 ) 


T + A 


3.2 x 10 3 




1.2 


Xbj + A 


1.7 x 10 3 (4.5 x 


10 2 ) 


6.2 x 10" 1 (1.6 x 10- 1 ) 



We conclude this section by assessing the theoretical uncertainties arising from the 
lack of knowledge of the precise values of the bottomonium MEs and from the freedom 
in the choice of the renormalization and factorization scales. Besides the QCD-motivated 
potential [30], the authors of Ref. [25] also employed three other phenomenological ansatze, 
namely, a power-law potential [33], a logarithmic potential [34], and a Coulomb-plus- 
linear potential (Cornell potential) [35]. Taking as the theoretical errors the maximum 
upward and downward deviations with respect to the results obtained with the QCD- 



motivated potential, we have /(9 T ( 1S ) 3 S[^ 



-0.93 

-o.oo 



9.28 



-10.85 
-2.70 



GeV 3 and /0*w>( ip ) [ 3 P n {1) 



GeV . Since there are no errors assigned to the CO MEs quoted in Ref. 



13], 



2.03 

for the purpose of this error analysis, we adopt the older results from Ref. [5], which are 

(see also the discussion in the 



3p(8) 



m 2 b 



T(1S) is* 



7 (8) 



fitted to earlier CDF data [36]. In compliance with Ref. [13 
context of Eq. (4.1) in Ref. [5]), we assume that ^0 T ( 15 ) 
We then repeat the evaluation of Table 2 with this set of bottomonium MEs, varying one 
ME at a time and combining the individual shifts in quadrature thereby allowing for the 
upper and lower half-errors to be different. The outcome is presented in Table 4. From 
Table 4, we observe that the NRQCD results for C = T(IS'), x&j(1P) carry errors of ±84% 
and ±29%, respectively, while the CSM result for C = XbjO-P) ma Y be increased by up 
to 47%. We note in passing that CO MEs of Ref. [13] fall outside the error bars of their 



10 



counterparts of Ref. [5], so that the NRQCD results given in Tables 2 and 4 are not actually 
mutually consistent. While the face values of the NRQCD results in Table 4 might by 
now be obsolete, we believe that the relative errors should still serve as a useful indicator. 
The CDF data [36] from which the CO MEs of Ref. [5] were extracted are only based 
on an integrated luminosity of 16.6 pb _1 collected in 1992-1993. The design value of the 
integrated luminosity to be delivered by the Tevatron until the end of 2009 is currently 
quoted as 8 fb _1 [37]. The combined data sample of the CDF and DO collaborations 
would then be three orders of magnitude larger than the one used in Ref. [36], leading 
to a reduction in statistical error by a factor of about 30. The experimental errors in 
the cross section measurements of bottomonium hadropro duct ion at the Tevatron should 
then be dominated by systematics. 



Table 4: Same as in Table 2, but using the CO MEs of Ref. [5] and taking into account 
the theoretical errors on the CS [25] and CO [5] MEs. 



C + D 


LHC 




Tevatron 




T + h 


(4.4 ±3.7) x 10 1 




(8.1 ±6.8) x 10- 1 




XbJ + h 


(2.1 ±0.6) x 10 3 (l.9+J;J x 10 2 ) 


(3i 


,±1.1) x 10 1 (3.6;! 


-1.7\ 
-0.0 J 


T + A 


(4.4 ±3.7) x 10 1 




(8.2 ±6.8) x 10- 1 




Xbj + A 


(2.1 ±0.6) x 10 3 (l.9+J;J x 10 2 ) 


(3i 


,±1.1) x 10 1 (3.6^ 


-1.6 \ 
-0.0 J 



Another plausible choice of renormalization and factorization scales is fi r = fif = mo- 
For pt <C ttlo {pt ^> thd), this scale choice is larger (smaller) than our default setting. 
The pt distributions are particularly sensitive to this change. In order to illustrate this, 
we recalculated the pr distributions shown in the upper panels of Figs. 1(a) and 2(a) using 
this alternative scale choice. The results are displayed in the upper and lower panels of 
Fig. 3, respectively. From Fig. 3, we observe that the px distributions are appreciably 
reduced in the lower px range, for p? ^ 25 GeV, as we pass from our default scale choice 
to the new one. This feature is more pronounced for the Tevatron than for the LHC. On 
the other hand, towards the upper end of the pr values considered in the case of the LHC, 
the shifts in cross section induced by switching to the new scale choice are insignificant. 
Such scale uncertainties are generally expected to be reduced as next-to-leading-order 
corrections are taken into account. 

4 Conclusions 

We studied the associated production of heavy quarkonia C and Higgs bosons D in photon- 
photon, photon-hadron, and hadron-hadron collisions to LO in the NRQCD factorization 
formalism. We considered all experimentally established heavy quarkonia, with 2S+1 Lj = 
XX 1 ^!, 3 ^/ (J = 0,1,2), and the SM Higgs boson D = H as well as the MSSM 
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Higgs bosons D = h,H, A, H ± . We listed the cross sections of all contributing partonic 
subprocesses, except for those with D = iJ ± , which are either suppressed by small Yukawa 
couplings or involve an additional heavy particle with a distinct signature in the final state. 
We presented numerical results for any combination of C — Y(IS'), xw(l-P) and D = h,A 
appropriate for the LHC and Run II at the Tevatron, concentrating on the region of the 
currently allowed MSSM parameter space characterized by a large value of tan (3 and small 
values of and uia- 

Criteria to discriminate NRQCD against the CSM include the following: (i) T + h 
and T + A associated production is allowed in NRQCD, but forbidden to LO in the CSM; 
and (ii) in NRQCD, the cross sections of Xbj + h and XbJ + A- associated production are 
a factor of 3 to 4 larger than in the CSM. 

In our MSSM benchmark scenario, there will be about 3000 (1) direct T(IS') + h — > 
l + l~ + bb events per year at the LHC (Tevatron) and similarly for the A boson. The 
Xbj — ► T(lS)j feed-down channels will add about 1500 (0.5) events to each of these 
yields. Results in the same ballpark are expected for the bottomonia with n = 2,3. We 
conclude that these signals should be clearly visible at the LHC, while they should provide 
a challenge for Run II at the Tevatron. 
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A Partonic cross sections 

In this appendix, we list the differential cross sections da/dt for processes (13)-(17) with 
D = h,A. The results for D — H,H may be obtained from those for D = h by accordingly 
replacing ijq and m h . Expressions of the partonic Mandelstam variables s, t, and u in 
terms of px, yc, an d Ud m ay be found in Eq. (4) of Ref. [21]. The mass m c — M of the 
heavy quarkonium is taken to be M = 2mg. 

da , — n mi x 2tt«2 f^) 2 M(2m? - s - t - uf 

dt {99 W[^0 \n) 9( m 2_ s _ i )2( m 2_ s _ u )2( 2m 2_t_ u )2' ^ 

^ faa - QO [3P«1 ^ = -^s{y h Q ) 2 

dt ^ \- J 1 ' 135M 3 s 2 (m 2 - s - t)\m\ - s - u)\2m\ -t-u)* 

xFj, J = 0,1, 2, (A.2) 
F = -10M 2 {196m 20 - 140mjf[8s + 7{t + u)} 
+ ml 6 [2756s 2 + 4844s(t + u) + 35(59t 2 + 134£m + 59m 2 )] 

- 2m^ 4 [1920s 3 + 5092s 2 (t + u) + s(4329t 2 + 9934£m + 4329m 2 ) 
+ 70(t + w)(17t 2 + 50£m + 17m 2 )] 

+ 2m^ 2 [1680s 4 + 6012s 3 (t + u) + 4s 2 (1908t 2 + 4409£m + 1908m 2 ) 

+ 7s(t + M)(589t 2 + 1790tM + 589m 2 ) 

+ 2(407t 4 + 2537t 3 M + 4402tV + 2537£m 3 + 407m 4 )] 

- 2m^°[968s 5 + 4462s 4 (t + u) + s 3 (7559t 2 + 17434£m + 7559m 2 ) 
+ 2s 2 (t + w)(2983t 2 + 9210tw + 2983m 2 ) 

+ s(2245t 4 + 14659t 3 M + 25776fV + 14659£m 3 + 2245m 4 ) 

+ (t + M)(335t 4 + 2874t 3 M + 5930* V + 2874£m 3 + 335m 4 )] 

+ m^[740s 6 + 4412s 5 (t + u) + 2s 4 (4803t 2 + 10882£m + 4803m 2 ) 

+ 2s 3 (t + M)(4959t 2 + 14882tM + 4959m 2 ) 

+ s 2 (5219t 4 + 34390t 3 M + 60782tV + 34390tM 3 + 5219m 4 ) 

+ Us{t + M)(99t 4 + 934t 3 M + 1994* V + 934tM 3 + 99m 4 ) 

+ 10(16t 6 + 239t 5 M + 1007tV + 1592tV + 1007tV + 239£m 5 + 16m 6 )] 

- 2m 6 [88s 7 + 728s 6 (t + u) + s 5 (2089t 2 + 4574tM + 2089m 2 ) 
+ 26s 4 (t + M)(107t 2 + 292tM + 107m 2 ) 

+ 5s 3 (375t 4 + 2270t 3 M + 3942tV + 2270£m 3 + 375m 4 ) 

+ s 2 (t + M)(643t 4 + 5906t 3 M + 12746* V + 5906*m 3 + 643m 4 ) 

+ s(lllt 6 + 1965t 5 M + 8844* V + 14224tV + 8844* V + 1965*m 5 + 111m 6 ) 

+ 10(* + m)(* 2 + 5*M + M 2 )(t 4 + 19* 3 M + 44*V + 19*M 3 + M 4 )] 

+ m^[20s 8 + 288s 7 (* + u) + 8s 6 (149* 2 + 314*m + 149m 2 ) 
+ 2s 5 (* + m)(1077* 2 + 2578*m + 1077m 2 ) 
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+ 2s 4 (969t 4 + 4952£ 3 m + 8206tV + 4952£m 3 + 969m 4 ) 

+ 2s 3 (t + w)(437t 4 + 3048t 3 M + 6210iV + 3048tM 3 + 437m 4 ) 

+ s 2 (179t 6 + 2630t 5 M + 1168HV + 18788tV + 11681tV + 2630tM 5 + 179m 6 ) 

+ 2s(t + u)(7t 6 + 258t 5 M + 1740* V + 3326tV + 1740* V + 258*m 5 + 7m 6 ) 

+ t 8 + 52* 7 m + 568tV + 2144* 5 m 3 + 3290*V + 2144tV + 568tV + 52*m 7 + m 8 ] 

- 2m 2 [12s 8 (t + u) + s 7 (91t 2 + 186tM + 91m 2 ) + 2s 6 (t + M)(121t 2 + 260tM + 121m 2 ) 
+ s 5 (308t 4 + 1361t 3 M + 2142tV + 136Hm 3 + 308m 4 ) 

+ 2s 4 (t + M)(101t 4 + 51H 3 m + 919tV + 511tM 3 + 101m 4 ) 

+ s 3 (65t 6 + 566t 5 M + 2070t 4 M 2 + 3182t 3 M 3 + 2070tV + 566£m 5 + 65m 6 ) 

+ s 2 {t + M)(8t 6 + 119t 5 M + 744tV + 1418tV + 744tV + 119tM 5 + 8m 6 ) 

+ stu(13t 6 + 204t 5 M + 852t 4 M 2 + 1334tV + 852tV + 204£m 5 + 13m 6 ) 

+ tu(t + M)(t 2 + 5tu + u 2 ){t A + 16t 3 M + 36tV + 16tM 3 + m 4 )] 

+ 9s 8 (t + m) 2 + 42s 7 (t + m) 3 + s 6 (79t 4 + 322t 3 M + 490iV + 322tM 3 + 79m 4 ) 

+ 4s 5 (t + m)(19£ 4 + 81t 3 M + 131tV + 81tM 3 + 19m 4 ) 

+ s 4 (39t 6 + 258t 5 M + 755t 4 M 2 + 1080t 3 M 3 + 755t 2 M 4 + 258£m 5 + 39m 6 ) 

+ 2s 3 (t + M)(5t 6 + 36t 5 M + 144* V + 246tV + 144tV + 36£m 5 + 5m 6 ) 

+ s 2 (t 8 + 10t 7 M + 91t 6 M 2 + 356tV + 552tV + 356tV + 91t 2 M 6 + 10tM 7 + m 8 ) 

+ 4stV(t + u)(t 2 + 5tM + M 2 )(3t 2 + 7tu + 3m 2 ) 

+ t 2 M 2 (t + M) 2 (t 2 + 5tM + M 2 ) 2 }, (A.3) 
Fx = -5(2m 2 -s-t- u) 2 {Amfs - 4m* 4 [10s 2 + 4s(t + «) + (*- u) 2 ] 
+ 2m^ 2 [62s 3 + 72s 2 (t + u) + s(23t 2 + lOfu + 23m 2 ) + 8{t + u)(t- u) 2 } 

- m^°[176s 4 + 368s 3 (t + u) + 2s 2 (119t 2 + 206£m + 119m 2 ) 

+ 2s(t + M)(47t 2 - 38tM + 47m 2 ) + (t - m) 2 (25£ 2 + 58£m + 25m 2 )] 
+ m 8 h [12As 5 + 400s 4 (t + u) + 4s 3 (109t 2 + 234£m + 109m 2 ) 

+ 4s 2 (t + M)(61t 2 + 68tM + 61m 2 ) + s(113t 4 + 86t 3 M - 118tV + 86tM 3 + 113m 4 ) 
+ (t + u)(t- M) 2 (19t 2 + 62tM + 19m 2 )] 

- m^Os 6 + 192s 5 (t + u) + 8s 4 (4K 2 + 96tM + 41m 2 ) 
+ 8s 3 (t + M)(34t 2 + 79tM + 34m 2 ) 

+ 2s 2 (83t 4 + 194t 3 M + 226* V + 194*m 3 + 83m 4 ) 
+ 4s(t + M)(19t 4 + 17t 3 M - 44t 2 M 2 + 17tM 3 + 19m 4 ) 
+ (t - M) 2 (7t 4 + 54t 3 M + 98t 2 M 2 + 54tM 3 + 7m 4 )] 
+ m£[4s 7 + 32s 6 {t + u) + 2s 5 (47t 2 + lUtu + 47m 2 ) 
+ 4s 4 (t + M)(31t 2 + 92tM + 31m 2 ) 

+ 2s 3 (51t 4 + 206t 3 M + 340tV + 206tM 3 + 51m 4 ) + 10s 2 (t + M) 3 (7t 2 - 2tu + 7m 2 ) 
+ s(27t 6 + 100t 5 M - t 4 M 2 - 140tV - t 2 M 4 + lOOtM 5 + 27m 6 ) 
+ (t + u)(t- u) 2 (t 4 + 16t 3 M + 38t 2 M 2 + 16tM 3 + m 4 )] 



14 



- m 2 h [6s 6 (t + uf + 2s 5 (t + u)(9t 2 + 26tu + 9m 2 ) 
+ s 4 (25t 4 + 124t 3 M + 2Ut 2 u 2 + 124£m 3 + 25m 4 ) 
+ 4s 3 (t + u)(6t A + 17t 3 u + 36tV + 17tu 3 + 6m 4 ) 

+ s 2 (15t 6 + 56t 5 M + 63tV + 60t 3 M 3 + 63tV + 56£m 5 + 15m 6 ) 

+ 2s(t + u)(2t 6 + Ut 5 u - 3tV - 18t 3 M 3 - 3tV + 14£m 5 + 2m 6 ) 

+ tu(t 2 - M 2 ) 2 (2t 2 + 9tu + 2m 2 )] 

+ s 5 (t 4 + 6t 3 M + 6fV + 6tM 3 + m 4 ) 

+ s 4 (t + u)(3t 4 + 8t 3 u + 18tV + 8tM 3 + 3m 4 ) 

+ s 3 (3t 6 + 8t 5 M + 25t 4 M 2 + 32tV + 25tV + 8£m 5 + 3m 6 ) 

+ s 2 (t + M)(t 6 + 6t 5 M + 3tV + 4tV + 3tV + 6tu 5 + m 6 ) 

+ stM(4t 6 + 5t 5 M - 2tV - lOtV - 2tV + 5tM 5 + 4m 6 ) 

+ t 2 M 2 (t-M) 2 (t + M) 3 }, (A.4) 

F 2 = 448mf - 16m 2 h ° [149s + 168(t + «)] 
+ 16m^ 8 [354s 2 + 800s(t + u) + 7(63t 2 + IMtu + 63m 2 )] 

- Amf [2111s 3 + 6834s 2 (t + u) + s(7491t 2 + 15766£m + 7491m 2 ) 
+ 140(i + M)(19t 2 + 46tM + 19m 2 )] 

+ 4m^ 4 [2358s 4 + 9272s 3 (t + u) + s 2 (14369t 2 + 29782tM + 14369m 2 ) 

+ s(t + M)(10041t 2 + 23470tM + 10041m 2 ) 

+ 7(363t 4 + 1796t 3 M + 2882tV + 1796tM 3 + 363m 4 )] 

- 2m^ 2 [4130s 5 + 18820s 4 (t + u) + 5s 3 (7093t 2 + 14430tM + 7093m 2 ) 
+ 26s 2 (t + M)(1329t 2 + 2954tM + 1329m 2 ) 

+ s(16995t 4 + 80520t 3 M + 127418tV + 80520£m 3 + 16995m 4 ) 

+ 56 (t + M)(57t 4 + 332t 3 M + 566tV + 332tM 3 + 57m 4 )] 

+ ml [5232s 6 + 28832s 5 (t + u) + 2s 4 (32427t 2 + 64870£m + 32427m 2 ) 

+ 2s 3 (t + M)(38653t 2 + 80838tM + 38653m 2 ) 

+ s 2 (52313t 4 + 233504t 3 M + 362542tV + 233504tM 3 + 52313m 4 ) 

+ 4s(t + M)(4710t 4 + 25333t 3 M + 41834* V + 25333tM 3 + 4710m 4 ) 

+ 7(t + M) 2 (377t 4 + 2748t 3 M + 5174tV + 2748tM 3 + 377m 4 )] 

- m 8 [2116s 7 + 15080s 6 (t + u) + 224s 5 (191t 2 + 377tM + 191m 2 ) 
+ 4s 4 (t + M)(15689t 2 + 31000tM + 15689m 2 ) 

+ s 3 (52571t 4 + 218938t 3 M + 332542tV + 218938£m 3 + 52571m 4 ) 

+ s 2 (t + M)(25783t 4 + 124384t 3 M + 197458tV + 124384tM 3 + 25783m 4 ) 

+ s(6841t 6 + 57484t 5 M + 171535t 4 M 2 + 241720t 3 M 3 

+ 171535tV + 57484tM 5 + 6841m 6 ) 

+ 7(t + M) 3 (99t 4 + 980t 3 M + 2162tV + 980tM 3 + 99m 4 )] 

+ m 6 [472s 8 + 4752s 7 (t + u) + 8s 6 (2279t 2 + 4436tM + 2279m 2 ) 
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+ 4s 5 (t + m)(8689£ 2 + 16336tw + 8689m 2 ) 

+ 2s 4 (18433t 4 + 72142t 3 M + 107398* V + 72U2tu 3 + 18433m 4 ) 

+ 4s 3 (t + m)(5706£ 4 + 24309£ 3 m + 37018* V + 24309tM 3 + 5706m 4 ) 

+ s 2 (8227t 6 + 60592t 5 M + 168433* V + 232088tV 

+ 168433tV + 60592*m 5 + 8227m 6 ) 

+ 2s(t + m)(785* 6 + 7956* 5 m + 25901tV + 37404tV 

+ 25901fV + 7956*m 5 + 785m 6 ) 

+ 35(f + m) 4 (3* 4 + 46t 3 M + 130tV + 46*m 3 + 3m 4 )] 

- m 4 h {Us 9 + 776s 8 (t + u) + 2s 7 (2229t 2 + 4262*m + 2229m 2 ) 
+ 8s 6 (t + m)(1453* 2 + 2593*m + 1453m 2 ) 

+ 4s 5 (4066t 4 + 15109t 3 M + 22095tV + 15109*m 3 + 4066m 4 ) 

+ 2s 4 (t + m)(6549* 4 + 25090t 3 w + 36846* V + 25090*m 3 + 6549m 4 ) 

+ s 3 (6151t 6 + 39452t 5 M + 101439*V + 136332*V 

+ 101439*V + 39452*m 5 + 6151m 6 ) 

+ s 2 (t + m)(1619* 6 + 13326* 5 m + 38633* 4 m 2 + 53884* 3 m 3 

+ 38633* 2 m 4 + 13326tM 5 + 1619m 6 ) 

+ s(207* 8 + 3214* 7 m + 16050* 6 m 2 + 38562tV + 50990* 4 m 4 

+ 38562*V + 16050t 2 M 6 + 3214*m 7 + 207m 8 ) 

+ 7(t + m) 5 (* 4 + 32t 3 M + 134*V + 32tM 3 + m 4 )] 

+ m 2 h [A8s 9 (t + m) + 2s 8 (271t 2 + 510*m + 271m 2 ) 

+ 2s 7 (* + m)(1023* 2 + 1726*m + 1023m 2 ) 

+ s 6 (3849* 4 + 13644* 3 m + 19606* V + 13644*m 3 + 3849m 4 ) 

+ 4s 5 (t + m)(1027* 4 + 3626* 3 m + 5160*V + 3626*m 3 + 1027m 4 ) 

+ 2s 4 (1290* 6 + 7467* 5 m + 17981* 4 m 2 + 23632* 3 m 3 + 17981* 2 m 4 + 7467*m 5 + 1290m 6 ) 

+ 10s 3 (* + m)(93* 6 + 626t 5 M + 1608*V + 2162tV + 1608*V + 626*m 5 + 93m 6 ) 

+ s 2 (173t 8 + 2052* 7 m + 8737* V + 19348* V + 24948*V 

+ 19348*V + 8737t 2 M 6 + 2052*m 7 + 173m 8 ) 

+ 4s(* + m)(3* 8 + 76* 7 m + 441* V + 1135t 5 M 3 + 1522*V 

+ 1135* 3 m 5 + 441*V + 76tM 7 + 3m 8 ) + 7tu(t + m) 6 (2* 2 + 17tu + 2m 2 )] 

- 24s 9 (* + m) 2 - 48s 8 (t + m)(3* 2 + 5*m + 3m 2 ) 

- s 7 (367* 4 + 1258t 3 M + 1786*V + 1258*m 3 + 367m 4 ) 

- s 6 (* + m)(515* 4 + 1708* 3 m + 2370*V + 1708*m 3 + 515m 4 ) 

- 2s 5 (215* 6 + 1162t 5 M + 2665* 4 m 2 + 3440* 3 m 3 + 2665* 2 m 4 + 1162*m 5 + 215m 6 ) 

- 2s 4 (t + m)(107* 6 + 630* 5 m + 1477* V + 1920* 3 m 3 + 1477* V + 630*m 5 + 107m 6 ) 

- s 3 (59* 8 + 558* 7 m + 2051*V + 4160* 5 m 3 + 5220*V 
+ 4160* 3 m 5 + 2051* 2 m 6 + 558*m 7 + 59m 8 ) 
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- s 2 {t + u)(7t 8 + 112t 7 u + 505tV + 1152tV + 1496tV 
+ 1152tV + 505tV + I12tu 7 + 7m 8 ) 

- stu(t + u) 2 (12t 6 + 85t 5 M + 242tV + 326£ 3 m 3 + 242tV + 85tu 5 + 12m 6 ) 

- 7t 2 u 2 (t + u) 7 , 

da / „— ri„rsM . \ 15 da / ^-^ 

99 -> <2Q 



_> gg [iffW] /, 



dt 

da 
~di 

x {9m* 2 



8 dt 



(A.5) 
(A.6) 



<? 5 - gg 3s} 



3 C (8) 



7rcr 



4Ms 3 (m^-s-t) 2 i 



2 (2 



it 



9m^ [2s + 3(f + «)] + m 8 [10s 2 + 40s(t + u) + 3(llt 2 + 23tM + 11m 2 )] 



- m 6 h [s 3 + 15s 2 (t + m) + 4s(9r 
+ m^[s 3 (t + u) + 2s 2 (5t 2 + 6tu 



■17tu + 9u 2 ) + 3(t + w)(7t 2 + 16fu + 7m 2 )] 
5m 2 ) + 4s(t + M)(4t 2 + 7tu + 4m 2 ) 



+ 7t 4 + 35t 6 u + 51tV + 35tM 3 + 7m 4 ] 

- m 2 h [s 3 (t 2 -tu + u 2 ) + s 2 (t + u)(3t 2 + tu + 3m 2 ) 

+ s(3t 4 + 14t 3 M + 16tV + 14tM 3 + 3m 4 ) + (t + u){t 2 + tu + u 2 ){t 2 + 7tu + u z 

+ tM[fi 2 (t 2 + M 2 ) + 2s(t + M)(t 2 + M 2 ) + {f + tU + M 2 ) 2 ]}, 



(A.7) 



da / 

dt ^ W 'J ^ = M 3 s 3 (m 2 - s - t) 3 (m 2 - s - M) 3 (2m 2 - t - m) 4 
x {36mf - Aml 8 [38s + 45(i + «)] 
+ m* 6 [300s 2 + 668s(t + u) + 3(13H 2 + 278£m + 131m 2 )] 



u 



s(639t 2 + 1330tM + 639m 2 ) 



473t 3 M + 750tV + 473tM 3 + 97m 4 )] 
3722tM + 1841m 2 ) 



-2m^ 4 [182s 3 + 586s 2 (t 
+ 6(t + M)(41t 2 + 98tM + 41m 2 )] 

+ m* 2 [308s 4 + 1256s 3 (t + u) + s 2 (1979t 2 + 4050£m + 1979m 2 ) 
+ 7s(t + M)(199t 2 + 450tM + 199m 2 ) + 4(97t 4 

- m^°[204s 5 + 940s 4 (t + u) + s 3 (1841t 2 
+ 2s 2 (t + M)(937t 2 + 2032tM + 937m 2 ) 

+ 7s(135t 4 + 622t 3 M + 966tV + 622£m 3 + 135m 4 ) 

+ 6(t + M)(33t 4 + 190t 3 M + 310tV + 190tM 3 + 33m 4 )] 

+ m|[104s 6 + 552s 5 (t + u) + 2s 4 (597t 2 + 1198tu + 597m 2 ) 

+ 4s 3 (t + M)(367t 2 + 776tM + 367m 2 ) 

+ s 2 (1074t 4 + 4779t 3 M + 7334tV + 4779tM 3 + 1074m 4 ) 

+ s(t + M)(407t 4 + 2145t 3 M + 3352* V + 2145tM 3 + 407m 4 ) 

+ 2(32t 6 + 303t 5 M + 915tV + 1280£ 3 m 3 + 915t 2 M 4 + 303£m 5 + 32m 6 )] 

- 2m 6 [16s 7 + 116s 6 (t + u) + s 5 (297t 2 + 602tM + 297m 2 ) 
+ 2s 4 (t + M)(197t 2 + 418tM + 197m 2 ) 

+ s 3 (335t 4 + 1481t 3 M + 2258t 2 M 2 + 148Hm 3 + 335m 4 ) 
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+ 2s 2 {t + m)(93£ 4 + 473£ 3 m + 723tV + 473tM 3 + 93m 4 ) 

+ s(54t 6 + 464t 5 M + 1319t 4 M 2 + 1814fV + 1319tV + 464tM 5 + 54m 6 ) 

+ 2(t + w)(3t 6 + 40t 5 M + 134t 4 M 2 + 186£ 3 m 3 + 134t 2 M 4 + 40tM 5 + 3m 6 )] 

+ m A h [As 8 + 52s 7 (t + u) + s 6 (189t 2 + 394tM + 189m 2 ) 

+ 3s 5 (t + w)(101t 2 + 230tw + 101m 2 ) 

+ 2s 4 (135t 4 + 626t 3 M + 960tV + 626£m 3 + 135m 4 ) 

+ s 3 (t + w)(165t 4 + 886t 3 M + 1330tV + 886tw 3 + 165m 4 ) 

+ s 2 (72t 6 + 623t 5 M + 1733tV + 2368tV + 1733* V + 623tM 5 + 72m 6 ) 

3 m 3 + 559tV + 185tM 5 + 16m 6 ) 
490tV + 376tV + 160tV + 28tu 7 + m s 



+ s(t + M)(16t 6 + 185t 5 M + 559tV + 768^ 



+ t 8 + 28t 7 u + 160tV + 376t"v 



5„.3 



-m 2 [4/(t 



M 



s 7 (27t 2 + 58tM 



2s 6 (t + M)(31t 2 + 80tM + 31m 5 



+ s 5 (67t 4 + 354t 3 M + 562tV + 354tM 3 + 67m 4 
+ 4s 4 (t + M)(10t 4 + 66t 3 t 



66tM 3 + 10m 4 



+ s 3 (t + M) 2 (17t 4 + 153t 3 M 
+ 2s 2 (t + M)(3t 6 + 42t 5 M + 
+ s(t + M) 2 (t 6 + 2H 5 m + 70tV 



- lOltV 

+ 213tV + 153tM 3 
122tV + 168t 3 M 3 4 



17m 4 



94t 3 M 3 



122tV + 42tM 5 + 3m 6 ) 



70tV + 21tM £ 



M 



+ 2tu(t + M)(t 2 + tU + M 2 )(t 4 + 8t 3 M + 12* V + 8tM 3 + m 4 )] 



+ (t + u)[s 8 (t + u)+ As 7 {f + 3tM + m 2 ) + 3s 6 (t + u)(2t 2 + 9tu + 2u 2 ) 
+ s 5 (4t 4 + 39t 3 M + 66* V + 39tM 3 + 4m 4 ) 
+ s 4 (t + m)(* 4 + 24t 3 M + 36* V + 24*m 3 + m 4 ) 

+ 2s 3 *m(* + m) 2 (6* 2 + 7tu + 6m 2 ) + s 2 tu(t + m)(5* 4 + 15* 3 m + 21* V + 15*m : 



5m 4 



stu(t 



111 



da / 



QQ 



(* 4 + 4* 3 M + 5*V + 4tM 3 + M 4 ) + *V(* + m)(* 2 + tU + M 2 ) 2 ]}, 

da 

J' 



3p 



(8) 
J 



*) = ?E(^ + Df(« 



QQ \ 3 PP] h) , 



(A.8) 
(A.9) 



da 
~dt 



x [2m 4 



(qq - QQ [ 3 S{ 8) 



y h Q 



27Ms 3 (2m 2 



t-uf 



2m 2 h {3s + t + u) + 4s 2 + 4s(t + u) + t 2 + m 2 



da 
~di 
x [2m 4 
da 
~dt 
da 
~dt 



da / 



(l9 
(l9 



QQ Pi 



1 D (8) 



h 



27M 3 s 3 (2m 2 h - t - u) 2 



2m 2 h (s + t + u) +t 2 + m 2 ], 



QQ [^P] h) 



16a da 



a. 



QQ 



3 D (8) 



h) 



dt 

2 



{gg - QQ h) , 



16a ln T/ , . da / 



J'=0 



QQ [ 3 P?\ h) , 



(A.10) 

(A.ll) 
(A.12) 

(A.13) 
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<Ja (77 - QQ [n] h) = ^%( 99 ^ QQ M h) , n = ^ , M 1} , (A. 14) 



V " ^ L J ' 9a 2 . df 

do-/ _n (1) i x 27ra2^^ 2 M 

dt ^ QQ i S ° \ A ) = 9s 2 (mi - s - t)\m\ - s - uf{2m\ - t - m) 2 

x {m^ - 2m A (s + t + u)+ m A [2s 2 + 2s(t + u) + t 2 + 4£m + m 2 ] - 2m 2 A tu(s + t + u) 

+ t 2 u 2 }, (A.15) 

^ fa, - OQ PP (1) 1 A) = -^(vq) 2 

dt ^ W ^ 3 J > \3hM*s\m\ - s - tY(m\ - s - u)\2m\ -t-u)* 

xGj, J = 0,1, 2, (A.16) 

G = -WM 2 {lQm A 6 (t - u) 2 - lQm 14 (t - u) 2 [3s + 4(t + u)] 

+ Am A 2 [s 4 - 4s 3 (t + u) + s 2 (17t 2 - 18tu + 17m 2 ) + 38s(t + u)(t - m) 2 

+ 2(t - u) 2 (13t 2 + 30tu + 13m 2 )] 

- 4m 1 4 °[4s 5 - 12s 4 (t + u) + 2s 3 (3t 2 - 2tu + 3m 2 ) + s 2 (t + u)(A7t 2 - A6tu + 47m 2 ) 
+ s(t - u) 2 (A7t 2 + UOtu + 47m 2 ) + 2(t + u)(t - m) 2 (11£ 2 + 34£m + 11m 2 )] 

+ m^[24s 6 - 48s 5 (t + u) - 72s 4 (t 2 + tu + m 2 ) + 20s 3 (t + M)(7t 2 + 2tu + 7m 2 ) 

+ 3s 2 (71t 4 + 88t 3 M + 2tV + 88tM 3 + 71m 4 ) + 6s(t + u)(t - u) 2 (l9t 2 + 62tM + 19m 2 ) 

+ (t - w) 2 (41t 4 + 276t 3 M + 486tV + 276tM 3 + 41m 4 )] 

- 2m A [8s 7 - 8s 6 (t + u) - 2s 5 (23t 2 + 26£m + 23m 2 ) + 2s 4 (t + M )(5t 2 + \2tu + 5m 2 ) 
+ s 3 (93t 4 + 192t 3 M + 230tV + 192£m 3 + 93m 4 ) 

+ s 2 (t + M)(65t 4 + 96t 3 M - 2fV + 96tM 3 + 65m 4 ) 

+ s(t - M) 2 (17t 4 + 128t 3 M + 230tV + 128tM 3 + 17m 4 ) 

+ (t + u)(t- u) 2 (5t 4 + 52t 3 M + llOtV + 52tM 3 + 5m 4 )] 

+ m 4 A [As 8 - 12s 6 (3t 2 + 4tM + 3m 2 ) - 4s 5 (t + m)(7£ 2 + 2tM + 7m 2 ) 

+ 2s 4 (37t 4 + 96t 3 M + 136tV + 96tM 3 + 37m 4 ) 

+ 2s 3 (t + M)(57t 4 + 98t 3 M + 130tV + 98tM 3 + 57m 4 ) 

+ s 2 (47t 6 + 164t 5 M + 197* V + 144iV + 197* V + 164tM 5 + 47m 6 ) 

+ 4s(t + u)(t- uf(t 4 + 15t 3 M + 34* V + 15tM 3 + m 4 ) 

+ (t 2 - M 2 ) 2 (t 4 + 22t 3 M + 66tV + 22tM 3 + m 4 )] 

+ 2m\{t + M)[2s 7 (t + u) + 4s 6 (t 2 + tu + m 2 ) - s 5 (t + M)(5t 2 + 8tu + 5m 2 ) 

- s 4 (19t 4 + 42t 3 M + 58tV + 42tM 3 + 19m 4 ) 

- s 3 (t + M)(17t 4 + 23t 3 M + 36t 2 M 2 + 23tM 3 + 17m 4 ) 

- s 2 (5t 6 + 15t 5 M + 2H 4 m 2 + 14t 3 M 3 + 21tV + 15tM 5 + 5m 6 ) 

- 3stu(t + U)(t- ufif + 4tM + M 2 ) - tM(t 2 - M 2 ) 2 (t 2 + 6tu + M 2 )] 
+ (t + M) 2 (S + t + U ) 2 [s A (t + M) 2 + 2s 3 (t + M)(t 2 + M 2 ) + S 2 (t 2 + M 2 ) 2 

+ tV(t-M) 2 ]}, (A.17) 



19 



Gi = 5{64mf - 4&m%[9s + 8(t + «)] 

+ 4m^ 8 [325s 2 + 584s(t + «) + 4(63t 2 + 134*u + 63m 2 )] 

- 4m^ 6 [566s 3 + 1563s 2 (t + u) + s(1381t 2 + 2890£m + 1381m 2 ) 
+ 20(t + u)(m 2 + 46tM + 19m 2 )] 

+ 4m^ 4 [619s 4 + 2376s 3 (t + u) + s 2 (3267t 2 + 6758£m + 3267m 2 ) 

+ 9s(t + m)(209£ 2 + 478tM + 209m 2 ) 

+ 363t 4 + 1796t 3 M + 2882tV + 1796£m 3 + 363m 4 ] 

- 2m^ 2 [864s 5 + 4422s 4 (t + u) + s 3 (8557t 2 + 17582£m + 8557m 2 ) 
+ 4s 2 (t + m)(1951£ 2 + 4298tw + 1951m 2 ) 

+ s(3269t 4 + 15024t 3 M + 23558fV + 15024tM 3 + 3269m 4 ) 

+ 8(t + M)(57t 4 + 332t 3 M + 566tV + 332tM 3 + 57m 4 )] 

+ m\°[748s 6 + 5072s 5 (t + u) + 2s 4 (663H 2 + 13590£m + 6631m 2 ) 

+ 2s 3 (t + M)(8637t 2 + 18662tM + 8637m 2 ) 

+ s 2 (11747t 4 + 51776t 3 M + 80050* V + 51776£m 3 + 11747m 4 ) 

+ 4s(f + M)(946t 4 + 4773t 3 M + 7714* V + 4773*m 3 + 946m 4 ) 

+ (t + M) 2 (377t 4 + 2748f 3 M + 5174tV + 2748£m 3 + 377m 4 )] 

- m\\lUs 7 + 1716s 6 (t + u) + 16s 5 (378t 2 + 775£m + 378m 2 ) 
+ 4s 4 (t + M)(2697t 2 + 5798tM + 2697m 2 ) 

+ s 3 (10665t 4 + 46174t 3 M + 70962* V + 46174tM 3 + 10665m 4 ) 

+ s 2 (t + M)(5765t 4 + 27152t 3 M + 42606* V + 27152*m 3 + 5765m 4 ) 

+ s(t + m) 2 (1471* 4 + 8390t 3 M + 14046tV + 8390*m 3 + 1471m 4 ) 

+ (t + m) 3 (99* 4 + 980t 3 M + 2162tV + 980*m 3 + 99m 4 )] 

+ m 6 A [20s 8 + 304s 7 (t + u) + 80s 6 (19t 2 + 39*m + 19m 2 ) 

+ 4s 5 (t + M)(927t 2 + 2002tM + 927m 2 ) 

+ 2s 4 (2549t 4 + 11026t 3 M + 16952tV + 11026tM 3 + 2549m 4 ) 

+ 8s 3 (t + M)(512t 4 + 2357t 3 M + 3666tV + 2357*m 3 + 512m 4 ) 

+ s 2 (t + M) 2 (1829t 4 + 9506t 3 M + 15006tV + 9506*m 3 + 1829m 4 ) 

+ 2s(t + m)(187* 6 + 1640t 5 M + 4899tV + 6900* 3 m 3 + 4899* V + 1640*m 5 + 187m 6 ) 

+ 5(t + M) 4 (3t 4 + 46t 3 M + 130*V + 46*m 3 + 3m 4 )] 

- m A [20s 8 (t + m) + 2s 7 (89* 2 + 182*m + 89m 2 ) + 4s 6 (t + m)(157* 2 + 342*m + 157m 2 ) 
+ 4s 5 (301* 4 + 1317t 3 M + 2036tV + 1317*m 3 + 301m 4 ) 

+ 2s 4 (* + m)(689* 4 + 3194t 3 M + 5012tV + 3194*m 3 + 689m 4 ) 

+ s 3 (941t 6 + 6676t 5 M + 18033* V + 24588* 3 m 3 + 18033* V + 6676*m 5 + 941m 6 ) 

+ s 2 (* + m)(353* 6 + 2874* 5 m + 8015*V + 11020* 3 m 3 + 8015*V + 2874*m 5 + 353m 6 ) 

+ s(57* 8 + 738* 7 m + 3246*V + 7294*V + 9442*V 

+ 7294fV + 3246*V + 738*m 7 + 57m 8 ) 
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+ (t + u) 5 (t A + 32t 3 M + 134tV + 32tu 3 + m 4 )] 

+ m 2 A [6s 8 (t + uf + 2s 7 (t + u)(19t 2 + 42tw + 19m 2 ) 

+ s 6 (lllt 4 + 500£ 3 m + 782tV + 500tM 3 + 111m 4 ) 

+ 4s 5 {t + w)(47t 4 + 225t 3 M + 364tV + 225tM 3 + 47m 4 ) 

+ 2s 4 (96t 6 + 691t 5 M + 1909tV + 2624tV + 1909tV + 691£m 5 + 96m 6 ) 

+ 2s 3 {t + m)(57£ 6 + 464t 5 M + 1294tV + 1778tV + 1294tV + 464tw 5 + 57m 6 ) 

+ s 2 (35t 8 + 440t 7 M + 181 It V + 3888tV + 4952tV 

+ 3888£ 3 m 5 + 181 It V 5 + 440£m 7 + 35m 8 ) 

+ 4s(t + M)(t 8 + 20t 7 M + 93tV + 215t 5 M 3 + 278tV 

+ 215tV + 93tV + 20tM 7 + m 8 ) + tu(t + M) 6 (2t 2 + I7tu + 2m 2 )] 

- (t + u)[s 7 (t + u)(t 2 + Atu + m 2 ) + s 6 (5t 4 + 28t 3 M + 50tV + 28tM 3 + 5m 4 ) 
+ 2s 5 (t + u)(5t A + 28t 3 u + 52tV + 28tw 3 + 5m 4 ) 

+ 2s 4 (t + M) 2 (5t 4 + 32t 3 M + 56tV + 32tM 3 + 5m 4 ) 

+ s s (t + M)(5t 6 + 56t 5 M + 160t 4 M 2 + 224t 3 M 3 + 160t 2 M 4 + 56£m 5 + 5m 6 ) 

+ s 2 (t 8 + 24t 7 M + 103tV + 224t 5 M 3 + 284tV + 224£ 3 m 5 + 103t 2 M 6 + 24£m 7 + m 8 ) 

+ stu(t + m)(4£ 6 + 19t 5 M + 46t 4 M 2 + 58t 3 M 3 + 46tV + 19£m 5 + 4m 6 ) 

+ t 2 M 2 (t + M) 6 ]}, (A.18) 

G 2 = 384m 2 4 2 - A8m^ [41s + 48(* + u)} 

+ 4m^ 8 [1071s 2 + 2604s(t + u) + 4(379t 2 + 802tM + 379m 2 )] 

- 4m^ 6 [1350s 3 + 4905s 2 (t + u) + s(6001t 2 + 12634£m + 6001m 2 ) 
+ 20(i + M)(115t 2 + 274tM + 115m 2 )] 

+ 4m^ 4 [1201s 4 + 5216s 3 (t + u) + s 2 (9677t 2 + 20274tM + 9677m 2 ) 

+ 6s(t + M)(1319t 2 + 3082tM + 1319m 2 ) 

+ 4(555t 4 + 2696t 3 M + 4298tV + 2696£m 3 + 555m 4 )] 

- 2m^ 2 [1768s 5 + 7602s 4 (t + u) + s 3 (17081t 2 + 35646£m + 17081m 2 ) 
+ 8s 2 (t + M)(2706t 2 + 6175tM + 2706m 2 ) 

+ 3s(4399t 4 + 20772t 3 M + 32874tV + 20772£m 3 + 4399m 4 ) 

+ 16(t + M)(177t 4 + 1000t 3 M + 1678tV + 1000£m 3 + 177m 4 )] 

+ m\° [2116s 6 + 8896s 5 (t + u) + 2s 4 (9939t 2 + 20686tM + 9939m 2 ) 

+ 2s 3 (t + M)(15713t 2 + 34582tM + 15713m 2 ) 

+ s 2 (30493t 4 + 138040t 3 M + 2 16542* V + 138040£m 3 + 30493m 4 ) 

+ Ss(t + M)(1812t 4 + 9605t 3 M + 15890tV + 9605£m 3 + 1812m 4 ) 

+ 2391t 6 + 21478t 5 M + 66169t 4 M 2 + 94100t 3 M 3 + 66169tV + 21478tM 5 + 2391m 6 ] 

- m^^s 7 + 4028s 6 (t + u) + 16s 5 (556t 2 + 1153tM + 556m 2 ) 
+ 4s 4 (t + M)(3631t 2 + 7486tM + 3631m 2 ) 

+ s 3 (18407t 4 + 76658t 3 M + 118270t 2 M 2 + 76658tM 3 + 18407m 4 ) 
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+ s 2 {t + m)(14251£ 4 + 67500t 3 M + 110090tV + 67500tM 3 + 14251m 4 ) 

+ s(5289t 6 + 43036t 5 M + 127559t 4 M 2 + 179784t 3 M 3 

+ 127559t 2 u 4 + 43036tw 5 + 5289m 6 ) 

+ {t + m)(645£ 6 + 7346t 5 M + 25323* V + 37052tV 

+ 25323tV + 7346tM 5 + 645m 6 )] 

+ m^lSSs 8 + 1104s 7 (t + u) + 8s 6 (357t 2 + 736tu + 357m 2 ) 

+ 8s 5 (t + m)(613£ 2 + 1150tw + 613m 2 ) 

+ s 4 (7042t 4 + 24892t 3 M + 36808* V + 24892tM 3 + 7042m 4 ) 

+ 8s 3 (t + M)(932t 4 + 3449t 3 M + 5458* V + 3449*m 3 + 932m 4 ) 

+ s 2 (4531t 6 + 30636t 5 M + 84881tV + 1 17904* V 

+ 84881*V + 30636*m 5 + 4531m 6 ) 

+ 6s(t + m)(209* 6 + 1962* 5 m + 6259*V + 9092* 3 m 3 + 6259* 2 m 4 + 1962*m 5 + 209m 6 ) 
+ (* + m) 2 (101* 6 + 1664* 5 m + 6795*V + 10240* 3 m 3 + 6795*V + 1664*m 5 + 101m 6 )] 

- m\[16s 9 + 124s 8 (* + m) + 2s 7 (229* 2 + 470*m + 229m 2 ) 
+ 4s 6 (* + m)(261* 2 + 446*m + 261m 2 ) 

+ 4s 5 (464* 4 + 1325* 3 m + 1794*V + 1325*m 3 + 464m 4 ) 

+ 2s 4 (* + m)(1259* 4 + 2936* 3 m + 4152*V + 2936*m 3 + 1259m 4 ) 

+ s 3 (2151* 6 + 10188* 5 m + 23655* 4 m 2 + 31340* 3 m 3 + 23655* 2 m 4 + 10188*m 5 + 2151m 6 ) 

+ s 2 (* + m)(967* 6 + 6318* 5 m + 17545*V + 25076*V 

+ 17545* 2 m 4 + 6318*m 5 + 967m 6 ) 

+ s(179* 8 + 2446* 7 m + 11502* 6 m 2 + 27386*V + 36254*V 
+ 27386*V + 11502* V + 2446*m 7 + 179m 8 ) 

+ (* + m) 3 (7* 6 + 226* 5 m + 1225* 4 m 2 + 1884*V + 1225*V + 226*m 5 + 7m 6 )] 

+ m 2 A [Us 8 {t + m) 2 + 2s 7 (* + m)(39* 2 + 58*m + 39m 2 ) 

+ s 6 (249* 4 + 580* 3 m + 674* V + 580*m 3 + 249m 4 ) 

+ 4s 5 {t + m)(127* 4 + 174* 3 m + 166*V + 174*m 3 + 127m 4 ) 

+ 2s 4 (306* 6 + 935* 5 m + 1375*V + 1480* 3 m 3 + 1375*V + 935*m 5 + 306m 6 ) 

+ 2s 3 (* + m)(201* 6 + 738* 5 m + 1426*V + 1854*V + 1426*V + 738*m 5 + 201m 6 ) 

+ s 2 (125* 8 + 1044* 7 m + 3733* 6 m 2 + 7996*V + 10328* 4 m 4 

+ 7996*V + 3733*V + 1044*m 7 + 125m 8 ) 

+ 4s(* + m)(3* 8 + 62* 7 m + 308* V + 784*V + 1058* 4 m 4 

+ 784*V + 308* 2 m 6 + 62*m 7 + 3m 8 ) 

+ tu(t + m) 4 (14* 4 + 135* 3 m + 206*V + 135*m 3 + 14m 4 )] 

- (* + m)[s 7 (* + m)(7* 2 - 4*m + 7m 2 ) + s 6 (35* 4 + 28* 3 m - 2*V + 28*m 3 + 35m 4 ) 
+ 2s 5 (* + m)(35* 4 + 12* 3 m - 8*V + 12*m 3 + 35m 4 ) 

+ 2s 4 (* + m) 2 (35* 4 + 8* 3 m + 6*V + 8*m 3 + 35m 4 ) 
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+ s 3 (t + m)(35£ 6 + 104t 5 M + 132t 4 M 2 + 176t 3 M 3 + 132£ 2 m 4 + 104tM 5 + 35m 6 ) 

+ s 2 (7t 8 + 64t 7 M + 197tV + 420tV + 548* V + 420t 3 M 5 + 197£ 2 m 6 + Q4tu 7 + 7m 8 ) 

+ 3stu(t + m)(4£ 6 + 19£ 5 m + 54tV + 74tV + 54tV + I9tu 5 + 4m 6 ) 

+ t 2 u 2 (t + m) 4 (7£ 2 + Wtu + 7m 2 )], (A.19) 



lo(8) 



(gg -> gg 

gg [ 3 si 8) 



d(7 

x {9m^ 2 -m\°[20s 
-m 6 4 [3s 3 + 23s 2 (t 
+ m\[3s 3 (t + u) 



A 



A 



15 (icr 
~8~ ~dt 



(gg -> gg 



(A.20) 



4Ms 3 (m^ - s - t) 2 (m\ - s - u) 2 (2m 2 A -t-u) 2 
- 27 (t + u)] + m 8 4 [14s 2 + 46s(t + u) + 3(llt 2 + 23tM + 11m 2 ) 
u) + 42s(t + m) 2 + 3(t + w)(7t 2 + 16fu + 7m 2 )] 



2s 2 (7t 2 + 13tu + 7m 2 ) + 2s(t + u)(9t 2 + 20tM + 9m 2 



+ 7t 4 + 35t 3 M 



51tV 



35tw 3 + 7m 4 



- mjjs^r + tu + m ) + s 2 (t + M)(3r + 7tM + 3m 2 ) 



+ s(3t 4 + 18t 3 M + 26tV + 18£m 3 + 3m 4 ) + (t + u)(t 2 + tu + u 2 )(t 2 + 7tu + u 2 )} 

+ tu[s(t + u) + t 2 + tu + u 2 ] 2 }, (A.21) 



da / 
-dt^ 99 
x {3m^ 6 [4s 2 
-m\ 4 [28s 3 - 
+ m\ 2 [28s 4 - 



gg M 8) 



A 



M*sZ{m 2 A -s- tf{m 2 A -s- uf{2m\ -t-u) 4 
As(t + u) + 3(t + u) 2 ] 



76s 2 (t -} 
108s 3 (t 



u) + s(73t 2 + 142tM + 73m 2 ) + 36(t + u) 
h«) + 2s 2 (97t 2 + 204tM + 97m 2 ) 



+ s(t + m)(161£ 2 + 330tM + 161m 2 ) + 12(t + M) 2 (5t 2 + lltu + 5m 2 )] 



-m A u [20s b + 76s 4 (t 



u 



4s 3 (43t 2 + 97tM + 43M 2 



+ 2s 2 (t + M)(127t 2 + 302tM + 127m 2 ) 

+ s(181t 4 + 783t 3 M + 1192tV + 783tM 3 + 181m 4 ) + 18(t + m) 3 (3£ 2 + 8tu + 3m 2 )] 



363tV 



229tM 3 + 48m 4 



+ m 8 A [8s 6 + 40s 5 (t + u) + s 4 (79t 2 + 190tM + 79m 2 ) 
+ 2s 3 (t + u)(73t 2 + 206tM + 73m 2 ) + 4s 2 (48t 4 + 229t 3 M 
+ s(t + M)(117t 4 + 560t 3 M + 866tV + 560£m 3 + 117m 4 ) 
+ 2(t + M) 2 (14t 4 + 79t 3 M + 129tV + 79tM 3 + 14m 4 )] 

- m A [8s 6 (t + m) + s 5 (23t 2 + 58tM + 23m 2 ) + 2s 4 (t + M)(19t 2 + 66tu + 19m 2 ) 

+ 2s 3 (38t 4 + 201t 3 M + 340tV + 201tM 3 + 38m 4 ) 

+ 2s 2 (t + M)(45t 4 + 226t 3 M + 384t 2 M 2 + 226tM 3 + 45m 4 ) 

+ s(45t 6 + 337t 5 M + 934tV + 1272t 3 M 3 + 934tV + 337tM 5 + 45m 6 ) 

+ 4(t + M) 3 (2t 4 + 16t 3 M + 27t 2 M 2 + 16tM 3 + 2m 4 )] 

+ m 4 A [2s 6 (t + m) 2 + s 5 (t + M)(5t 2 + Utu + 5m 2 ) 

+ 2s 4 (7t 4 + 31t 3 M + 58t 2 M 2 + 31tM 3 + 7m 4 ) 
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+ s 3 (t + w)(29t 4 + 122t 3 M + 262t 2 u 2 + I22tu 3 + 29m 4 ) 

+ s 2 (27t 6 + 180t 5 M + 539t 4 w 2 + 768t 3 w 3 + 539tV + 180tw 5 + 27m 6 ) 

+ s(t + w)(10t 6 + 86t 5 u + 255tV + 342tV + 255tV + 86tw 5 + 10m 6 ) 

+ (t + uf{f + m b u + 63tV + 88tV + 63tV + 18tM 5 + m 6 )] 

- m 2 A [s 5 (t 4 + t 3 u + tu 3 + m 4 ) + 4s 4 (t + u)(t 4 + 2t 3 u + 5fV + 2tw 3 + u 4 ) 

+ 2s 3 (3t 6 + 16t 5 w + 48tV + 76tV + 48tV + 16tw 5 + 3m 6 ) 

+ 2s 2 (t + u)(2t 6 + 15t 5 w + 43tV + 72tV + 43tV + 15tw 5 + 2m 6 ) 

+ s(t 8 + 15t 7 w + 64tV + 142tV + 180tV + 142tV + 64tV + 15tw 7 + u 8 ) 

+ 2tu(t + u) 3 {t 4 + 5t 3 M + 6tV + 5tw 3 + m 4 )] 

+ tu(t + u)[s 4 {t + u){t 2 + tu + u 2 ) + s 3 (3t 4 + 6t 3 w + 14tV + 6£m 3 + 3m 4 ) 
+ s 2 (t + u)(3t 4 + At 3 u + 14* V + 4tM 3 + 3m 4 ) 

+ s(t + u) 2 {t 4 + 2t 3 M + 4iV + 2tM 3 + m 4 ) + tu(t + u){t 2 + tu + u 2 ) 2 }}, 



da 
~dt 



(gg - QQ [ 3 Pj 8) ] a) 
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da 



dt 

x [2m 4 A - 2m 2 A (s + t + u)+t 2 + it 2 ], 



8 B 2j,+ y 

27Ms 3 (2m 2 A -t- U y 



gg -> QQ 



3p(l) 



A 



(A.22) 
(A.23) 



(A.24) 



da 



a i{ q - q ->QQ[ 1 Pi 8) ]A) = 



-lQna 2 s (y$)'' 



27M 3 s 3 ( y 2m 2 A -t- U y 



x {8m 6 s - 2m 4 A {2s + t + u) 2 - 2m 2 A [s{t - m) 2 - (t + u) 3 ] - (t + u) 2 {t 2 + m 2 )}, (A.25) 



da / 



dt 
da 
~dl 
da 



Q Q[^}A) = ^^( 99 
1 J ' a s at v 



QQ 



A 



da / 



3p(8) 



dt 



(77 



QQ 
QQ [n] A) 



16a 



a. 



ls J'=0 

512a 2 da 



QQ 



3p(l) 



A), 



9a? rft 



(A.26) 
(A.27) 
(A.28) 
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pp -> C D + X at the LHC, C = T, x, 



P T [GeV] p T [GeV] 

10 20 30 40 50 60 10 20 30 40 50 60 




m CD [GeV] m CD [GeV] 

(a) 

Figure 1: (a) px distributions da/dpr (upper panel) and men distributions da/dmcD 
(lower panel) in fb/GeV, (b) yc distributions da /dye (upper panel) and y& distributions 
da/dyo (lower panel) in fb, and (c) total cross sections a in fb as functions of tan/3 
(upper panel) and m D (lower panel) of pp — > CD + X, where C = T(1S), Xw(l-P) and 
D = h (left columns) or D = A (right columns), at the LHC. The default values of 
the MSSM input parameters are tan/3 = 50 and mo = 100 GeV. It is summed over 
C = Xbo(l-P), Xm(I-P), Xb2^-P)- In each frame, the NRQCD (solid line) prediction for 
C = T(1S) and the NRQCD (dashed line) and CSM (dotted line) ones for C = Xbj(lP) 
are shown separately. 
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pp -> C D + X at the LHC, C = T, %, 




pp -> CD + X at the LHC, C = T, % 



tan(3 tan(3 
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(c) 

Fig. 1 (continued). 
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pp -> C D + X at the Tevatron, C = T, x bJ 
P T [GeV] p T [GeV] 
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Figure 2: Same as in Figs. l(a)-(c), but for pp — > CD + X in Run II at the Tevatron. 
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pp -> C D + X at the Tevatron, C = T, %, 




pp -> C D + X at the Tevatron, C = T, x, 
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Figure 3: Same as in the upper panels of Figs. 1(a) and Figs. 2(a), but for \x T = /if = m D . 
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